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Abstract 

We derive for the Navier-Stokes equation an exact equation satis- 
fied by the dissipation rate correlation function, {e{x + f,t + T)e{x,t)). 
In the equal time limit, for the homogeneous, isotropic state of fully- 
developed turbulence, we show that the correlation function behaves 
as Ar~^^ + Br~^^ with = 2 — C,q and ^2 = z'l — Q/^ for r in the inertial 
range; the ^'s are exponents of velocity structure functions and z'l is a 
dynamical exponent characterizing the 4th order dynamical structure 
function. This provides the first direct derivation of the exponents of 
the dissipation-rate correlation. 
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The statistical properties of the energy dissipation rate e{x, t) defined by 

t) = 7;Y.i + djUi f (1) 

where di = d/dxi have played a crucial role in our understanding of fully- 
developed turbulence in incompressible fluids, [jl], Q In the original Kolmogorov 
theory e is replaced by (e) and the spatial fluctuations are ignored; the effect 
of fluctuations of e pointed out by Landau have been explored in the context 
of the lognormal model and its multifractal generalizations. The intermit- 
tent behavior of turbulent fluctuations is reflected in the power-law behavior 
of the correlations of e: 

(e(f)e(x + f)) ~ {r/Ly (2) 

where L is a length scale characteristic of the large-scale flow and r = |r| 
belongs to the inertial range. Simple dimensional analysis, noting that the 
dimension of e is V^/L, yields the identiflcation fi = 2 — (q, the exponents 
of the g*^-order (longitudinal) structure function, (g are deflned by 

S, = {[5u-f]^) {r/Lf^ (3) 

where 5u = u{x + r,t) — u{x,t). Within the original Kolmogorov theory 
(q = 2 and consequently fi = 0. Thus the deviation of fi from zero is a mea- 
sure of the degree of intermittency and is an important quantity for under- 
standing fully developed turbulence. This breakdown of simple Kolmogorov 
scaling has been studied experimentally and a review of the experiments 
gives a "best" estimate for /x of 0.25 ± 0.05 which is consistent with the ex- 
perimentally measured value of (q. In this Letter we provide a simple and 
direct derivation of the values of /i for the Navier-Stokes equation by deriving 
the exact equation satisfled by the dissipation rate correlations. The equa- 
tion for the dissipation rate contains two contributions, one a second spatial 
derivative of an appropriate sixth-order structure function and another a 
second temporal derivative of a fourth-order structure function; there are, in 
addition, pressure-dependent terms and no other velocity-dependent terms. 
The use of dynamical structure functions, i.e., Sg in Eqn. deflned with 
Su = u{x + r,t + t) — u{x,t), is key to our derivation. Both spatial and 
temporal derivatives of the dynamic structure functions occur naturally and 
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the equal time limit of the derivatives is related to the correlation of the 
energy dissipation rates. We illustrate the method first with an application 
to the Id stochastic Burgers equation which provides an extremely fruitful 
model for elucidating some of the conceptual and mathematical features of 
3d turbulence. After providing the derivation in the Navier-Stokes case we 
comment on our results and discuss some experimental implications. 

Derivation for the stochastic Burgers equation: We discuss first the 
Burgers equation (describing a one-dimensional, compressible fluid without 
pressure) driven by a stochastic force f{x,t): 

du/dt + udu/dx = vd'^ujdx + / . 

The stochastic force is a Gaussian white-noise with algebraic spatial correla- 
tions given in fc-space by 

{f{k,t)f{k',t')) = Do\kf6k+k',oS{t-t') (4) 

with —1 < j3 < 0. For these values of (3, intermittent behavior of the 
structure functions has been shown to occur. We will first illustrate our 
approach by doing the simplest calculation. We employ the notation u = 
u{x, t), u' = u{x', t') where x = R + r/2, x' = R — r/2, t = T + r/2, and 
t' = T — r/2. We wish to find an equation for (ee') where e = v{du/dxY is 
the dissipation rate. Multiplying the Burgers equation by u we have 

vud'^u/dx + fu = udu/dt + u^du/dx 

and the corresponding equation for u' . The idea is simply to multiply the 
two equations to obtain 

[vud'^u/dx + fu] [vu'd'^u /dx' + fu] = [udu/dt + u^du/dx] [u'du /dt' + u'^du' /dx] 

(5) 

and average over the homogeneous, steady state of Burgers turbulence; the 
identity vud'^u/ dx"^ = vd'^{u'^ /2)/ dx"^ — e is crucial. We obtain after some 
straightforward rearrangements of the terms on the right-hand side using 
d/dt = d/d-r and d/dt' = —d/dr when acting on averages in the steady 
state. 
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All the terms in the above equation are assumed to be evaluated in the 
limit r — > O"*". The only reason for the approximate sign is that we have not 
displayed terms which depend explicitly on u or the noise and are negligible in 
the inertial range as z/ — 0: for example, a term such as u"^ 8^84 is negligible 
since is finite. In the equal-time limit the noise-velocity correlations that 
occur on the left-hand side of Eqn. can be evaluated by using the Donsker- 
Novikov-Varadhan result 0; of these the only term that does not vanish in 
the inertial range is the subtracted term (e)^ in Eqn. (^). 

We immediately observe that we obtain the second temporal derivative 
of a fourth-order structure function and the second spatial derivative of a 
sixth-order structure function which yield the two exponents referred to in 
the abstract. Of course, the equation is not manifestly form-invariant under 
Galilean transformations and we have to understand the role of the cross 
spatio-temporal derivative. We address this issue using a different version of 
the above equation. 

A more elegant form can be obtained by multiplying Burgers equations 
for u and u' hy 5u = u — u' = u{x, t) —u{x', t'), multiplying the two equations 
and averaging over the homogeneous turbulent state as before. We obtain 

Dv Dii' 

{ 6ui ud\ + f )6ui ud'\' + /') ) = {5u — 6u^) (7) 

where D/Dt represents the convective derivative, i.e., Du/Dt = du/dt + 
udu/dx, and we have used the shorthand notation d = d/dx and d' = d/dx'. 
Evaluating the purely viscous terms on the left-hand side of Eq. (^) yields 

{i6ufd\d'\') = -2(ee') + '^dtUdu)^) + z/ (e + e')(H' ) 

where the only non- vanishing term in the inertial range is — 2(ee') in the 
limit z/ — i> since as before the rest are products of powers of u and finite 
correlation functions. The noise terms can again be evaluated in the r ^ O"*" 
limit to yield 2(e)^ — {l/12)[dS^/ drY apart from terms which vanish as z/ — *■ 
0. The right-hand side of Eq. contains four terms and their evaluation is 
facilitated by kinematic results which can be obtained in a straightforward 
way such as 

dl{{5uf) = 12 {{5ufdtudt>u') (8a) 
dl{{5uf) = 30{{5uYd^ud,m') (86) 
d^{{u + u'f{6uy) = {[-2{6uy + 12{u + u'f{6uf]d,ud,>u') (8c) 
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These enable one to simplify the temporal and spatial second derivatives. 
The term involving cross (space and time) derivatives can be simplified by 
employing 

drdr { {6u)\ ) = ( [12{6u)\ + A{6u)^]d^u d^u ) (9a) 
drdr{{5uYu') = {[12{6u)\' - A{6uf]d^>u' dtu) (96) 
drdr{{6u)^) = 20{{6uf]d^udt'u) = 20{{6uf]d^m' dtu) . (9c) 

We obtain after a few manipulations 

(10) 

where we have used the notation 5*^,^ = {{u — u'Y {u + u'Y). Note that 
d'^ 5*4 + drdr Si^i + \d'^Si^2 is the convective second derivative of the fourth- 
order structure function, D'^S^^/ Dt"^. This term is manifestly Galilean-invariant 
as is evident if we recall that dr{ {u — u'Y) + {l/2)dr{ {u + u'){u — u'Y ) is 
Galilean- invariant. The last term on the right-hand side which arises from 
the noise- velocity correlations leads to r^''^"^ which yields for the dissipation- 
rate correlation exponent precisely the scaling value of 2 — C^] the explicit 
value is given by 2 — 2(^3 = 2 + 2(3 where the last equality follows from 
the von-Karman-Howarth relation derived previously by us[Q, C3 = — /5. In 
the multifractal regime this term is, of course, subdominant. The crucial 



feature of Eq. ( |10|) is that the non-scaling behavior of the dissipation rate 
correlations is determined by (1) d^S^ term which leads to /ii = 2 — (^g and 
{2)D'^Si/ Dt'^ . Let us define z'l to be the dynamical scaling exponent describ- 
ing the r ^ limit of the Galilean- invariant, convective second derivative 
of 5*4. This yields another possible inertial range behavior for (ee') with 
/i2 = z'l — C4. The need for a sequence of dynamical exponents is a conse- 
quence of the occurrence of temporal multiscaling in the dynamical structure 
functions as has been emphasized earlier 0; thus different order temporal 
derivatives of Sp{r, r) can lead to different dynamical exponents. In dealing 
with dynamic structure functions it is important to recall that we have used 
the Eulerian description; therefore, ordinary dynamic scaling and a fortiori, 
dynamic multifractality, are complicated by the presence of sweeping terms. 
For example, in Sp{r = 0, r) (obtained from measurements of velocity differ- 
ences at a given point at finite values of the time difference) the kinematic 
exponent z = 1 arising from sweeping occurs. [110] However, in the (Galilean 
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invariant) convective derivative which occurs in the equation above in the 
T ^ hmit only the intrinsic, dynamical exponent occurs. Thus the exact 
equation neatly picks out the intrinsic dynamical exponent. We will consider 
the implications of the two intermittency exponents after the corresponding 
derivation for the three-dimensional problem. 

Derivation for the Navier-Stokes equation: We consider the Navier- 
Stokes equation for the velocity field ^(a;, t) driven by a stochastic driving 
force f{x, t) with zero mean and variance given by 

( mW) ) = p^Ak) ^(^)%+fe^o - (11) 

where Pijik) is the transverse projection operator given by 5ij — {kikj/k'^). 
The noise covariance D(k) is assumed to be peaked around ko ~ 1/L with 
a narrow width. In contrast to the Burgers equation the detailed form of 
the noise correlation is not crucial in the 3d problem; the noise maintains 
a fully-developed turbulent state and allows one to define averages as noise 
ensemble averages. We will find it useful to define the quantity 

eij = udeUidiUj . (12) 

The dissipation rate e (cf. Eq. (1)) of an incompressible fiuid obeys the 
relation 

€ = eu- z/V'p (13) 

where p = p/p {p is the constant density) and the summation convention 
of summing over repeated indices is used. We remark that (ejj) oc 6ij in 
isotropic turbulence. 

We use the notation x = R + (l/2)f, t = T + (1/2)t and x' = 
R — (l/2)f, t' = T — (l/2)r. Multiplying the Navier-Stokes equation for 
Ui = Ui{x, t) by Ui and summing over i and using vUiV'^Ui = uV'^{v? /2) — en 
one finds 

- en + vV'^{v?/2) + f - u ^ UidtUi + UiUidiUi + Uidip . (14) 

Again we write a similar equation for u[ — Ui{x',t') and multiply the two 
equations and average over the homogenous, steady state of isotropic turbu- 
lence leading to 

(ene' -) - - ( Ciiti'^ e-jM^ ) — V^V^(mV^) noise terms 
2 4 
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1 (9^ 1 

I 2 /2 \ J-ty ,/\2/2\ -LL/ , 2 >2 \ 

+ pressure terms . (15) 

The second and third terms on the left-hand side are neghgible in the inertial 
range and will be suppressed hereafter; the velocity terms clearly have the 
same form as in the Burgers equation consisting of a second derivative with 
respect to r of the sixth-order structure function and a second derivative with 
respect to r of a fourth-order structure function. We have not displayed 
the noise and pressure terms since this calculation merely illustrates the 
simple steps involved in the derivation. We will discuss them in the Galilean- 
invariant form given next. 

The manifestly Galilean-invariant form is somewhat more difficult to ob- 
tain in the 3 — d problem. We proceed from the Navier-Stokes equation in two 
different ways and add the results appropriately. First multiply the Navier- 
Stokes equation for Ui and u[ by 6ui = Ui{x,t) — Ui{x',t') and multiply the 
two equations to obtain 

Dv ■ Du'- 
5u, {vV\, + /,) 5u, {vV% + j;) = 6u,{^ + diP) 6u, ( + d'^p' ) . 

(16) 

We now average this equation over the isotropic, homogeneous, steady state 
of turbulence. On the left-hand side we obtain, in addition to the noise terms, 
— ( ^ue'jj + eije'j^j ) as the only terms which survive in the inertial range. We 
do a similar set of manipulations with the Navier-Stokes equation for Ui 
multiplied by 6uj and obtain 

Dti ■ Dv' 

(17) 

We add Eqn. ([T6|) to twice Eqn. (p!7D and average the resulting equation over 
the turbulent state. The right-hand side of the resultant equation includes, 
apart from the pressure terms, 

Dui^ Du' Du. Du[. 

(Sui—^Sui—^) + 2(6ui—^6u.:—-^) . 
^ Dt ^ Dt' ' ^ ^ Dt ^ Dt'' 

We can show after some algebraic manipulations involving kinematic rela- 
tions which are the three-dimensional generalizations of those displayed in 



7 



Equations (8) and (9) that these terms correspond to 

{l/l&){d''/dridrj){5ui5uj{5u-5uf) - {l/A){Dy Dt^){{Su ■ SuY) 

where D^Dr^if) = / dr^ {f) + / dTdri{{u^+u[)f) + {l/A)d^ / dridrj{{u^+ 
u'j){uj + u'j)f). This shows the occurrence of the two key structure function 
derivatives. We have evaluated the pressure terms and find 

1 92 I JJ2 

2(eiie'- - + 2eye.) + noise terms = ——— {6ui6uj{6u ■ 6uf ) - --—{{5u-5uf) 

+ ( 25ui5ujdipd'S' + 5ui5uidjpd'3' — — — ( {dip + d[p')5ui5u ■ 5u) 

+ ^dr, ( SuiSu ■ 5udui{dip - d[p') ) . (18) 

We emphasize that apart from terms neghgible in the hmit z/ ^ the 
above equation is exact (the neglected terms have been evaluated.) To com- 
plete our discussion we must consider the noise and pressure terms in Eq. 
(PD. We first note that 

(Q,e^) = (ee') + uVl{tP' + t'p) + v^VlVl{PP') 

which is equal to (ee') in the inertial range. We expect the diagonal elements 
of (diUidjUi) to yield the most singular terms and these survive in the z/ — >• 
limit which cuts off the short-distance singularities. With this observation 
we see that only the terms with i = j contribute in {eije'^j); using the isotropy 
of the turbulent state we find that the left-hand side of Eq. ([T8|) reduces to 
(10/3) (ee') in the inertial range. The noise terms can be evaluated and the 
non-trivial term in the t —>■ limit is precisely (10/3) (e)^. Thus the left- 
hand side of Eq. (p!8| ) is (10/3) [(ee') — (e)^] in the inertial range. Finally, 
we observe that the pressure terms depend upon derivatives of the pressure 
dip. The kernel in the inversion of V^p = —diUjdjUi is Coulombic and 
long-ranged; however, the pressure-derivative terms can be written with a 
kernel which is dipolar and hence, shorter-ranged leading to more convergent 
integrals. Thus writing the pressure contributions in terms of velocities one 
concludes plausibly that these contributions will not be more dominant than 
those due to the terms D'^Sa/Dt'^ and d'^Se/dr'^. 

We discuss our results next. The behavior of (ee') is thus determined by 
the two terms in the limit r = 0, analogous to those in the Burgers equation. 
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The term dr^dr {Sui6uj{6u ■ 5uY) yields /ii = 2 — C^q. The dynamical term, 
D"^ /Dt'^{5u- Su)^), yields the identification /i2 = z'^ — (4. This provides a 
transparent derivation directly from the Navier- Stokes equation of the two 
dominant exponents characterizing dissipation-rate correlations, one which 
depends purely on the static structure function exponent (2 — (q) and the 
other which involves dynamical behavior {z'l — (4). 

We proceed further, motivated by a similar strategy in the theory of phase 
transitions [0, and make an ansatz that the two terms are equally dominant. 
This leads to the identification [|l^ 

2 - Ce = 4' - C4 . (19) 

This relation connects multifractality in spatial correlations to multifractality 
in temporal correlations. [^] We then obtain the Kolmogorov result for the 
leading dissipation rate correlation exponent, 2 — (q. 

An earlier paper by us|ll[] provided a justification for the two exponents 
for the dissipation rate correlations using detailed equations for various struc- 
ture functions explicitly, identifying terms in different equations and invoking 
plausible comparisons with the Burgers equation result. We note that other 



relations such as /i = 2^2 — C4 have been proposed in the literature. |jT^ Other 



derivations of the Kolmogorov relation include one based on fusion rules for 



equal time multipoint correlation functions by L'vov and Procaccia|]T5[ who 
also pointed out another scenario which yields fi = 2(2 — C4 as in ref. |]T3 . 

We point out that in our equation the subtracted correlation, (ee') — (e)^, 
arises naturally. In the inertial range we expect the term (L/r)^ to dominate 
over the constant (e)^ term. However, this might require r « L since ix 
is small and this renders the experimental determination of /z more difficult. 
However, for r « L whether one subtracts or not one should obtain the 
correct exponent as expected from these theoretical considerations, and this 
is clear in Ref. 0. 
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